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ABSTRACT 

Runaway stars ejected from the Galactic disk populate the halo of the Milky 
Way. To predict the spatial and kinematic properties of runaways, we inject 
stars into a Galactic potential, compute their trajectories through the Galaxy, 
and derive simulated catalogs for comparison with observations. Runaways have 
a flattened spatial distribution, with higher velocity stars at Galactic latitudes 
less than 30°. Due to their shorter stellar lifetimes, massive runaway stars are 
more concentrated towards the disk than low mass runaways. Bound (unbound) 
runaways that reach the halo probably originate from distances of 6-12 kpc (10- 
15 kpc) from the Galactic center, close to the estimated origin of the unbound 
runaway star HD 271791. Because runaways are brighter and have smaller ve- 
locities than hypervelocity stars (HVSs), radial velocity surveys are unlikely to 
confuse runaway stars with HVSs. We estimate that at most 1 runaway star 
contaminates the current sample. We place an upper limit of 2% on the fraction 
of A-type main sequence stars ejected as runaways. 
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Subject headings: Galaxy: kinematics and dynamics — Galaxy: structure — 
Galaxy: halo — Galaxy: stellar content — stars: early-type 



Introduction 



Radial velocity surveys of the Milky Way halo can now reveal rare velocity outliers 
resulting from dynamical processes throughout the Galaxy. Presently there are three known 
categories of v elocity outliers among main sequence stars: (1) runaway stars originating in the 
Galactic disk (IHumason fc Zwickvlll947l). ( 2) velocity outliers from tidal disruption of dwarf 
satellites within the h alo (Ibata et al. Il994 ). and (3) hypervelocity stars (HVSs) originating 
in the Galactic center ( Brown et al. 2005h. Recent simu lations predict the observable spatial 
and velocity distrib utions for HVSs (IKenyon et al.ll2008l ) and for remnants of tidal disruption 
(lAbadi et al.ll2008l ). Although runaway stars have been known for a long time, there has 
been only limited investigation of their observable spatial/velocity distribution throughout 
the Galaxy. 

Here, we simulate the production of runaways and trace their orbits through the Galaxy, 
concentrating on runaways that reach the halo. We derive the spatial and velocity distribu- 
tions of runaways as a function of mass and use these results to predict observable quantities 
for comparison with data from radial velocity surveys. From these results, we estimate the 
fraction of runaways in surveys of halo stars. We also evaluate the likelihood of runaway 
star contaminants in targeted surveys for HVSs. 



1.1. History and Challenges 



First reported by lHumason &: ZwickyI (119471 ). runaway stars are short-liv ed stars at un- 



expectedly large distari ces and velocities relative to their probable site of origin (IBlaauw fc Morgan 



(Gies 


1987; 


Stone 


1991) 



velocity dispersion, 30 km s~^, roughly 3 times larger than the velocity dispersion of non- 
runaway stars in the Galactic disk. Howeve r, several runaways have velocities ^ 100 km s~^ 
relative to the local standard of rest (e.g. 



Gies 


1987; 


Stone 


1991; 


Martin 


2006 


) and large 


(Greenstein & Sareent 


1974; 


Martin 


2006). The 



large velocity dispersion and vertical scale height suggests a dynamical process that ejects 
runaways from the thin disk into the halo. 



Runaways probably originate within star-forming regions in the disk of the Milky Way. 
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Disruption of massive binaries is the likely source of runaways. In the binary supernova 
mechanism, a main sequence runaway is ejected when its former companion star explodes 
as a supernova; subsequent mass loss an d the kick velocity from th e asymmetric explosion 
are sufficient to unbind the binary (see iBlaauwl Il96ll : iHilld Il983l ). In this scenario, the 



runaway moves away with roughly the sum of the kick velocity and the orbital velocity of 
the binary, a velocity physically limited by the orbital velocity at the surface of the stars. 
In the dynamical ejection mechanism, runaways are ejected in dynamical t hree- or four- 



body interactions; the outcome is any co mbination of single stars and binaries (jPoveda et al. 



19671 : lHofferlll983l : iHut fc Bahcalll Il983l ). The maximum ejection velocity from d ynamical 
binar y-binary encounters is formally the escape velocity of the most massive star (ILeonard 
199lh . 



Understanding whether the proposed formation mechanisms can produce the observed 
runaway star populations requires high quality kinematic data and a clear understanding of 
the dynamics of close binaries and dense star clusters. For large ensembles o f runaways, the 



low b inary frequency of runaways favors the dynamical ejection mechanism ( Gies fc Bolton 



19861) . However, the pred icted kick velocity of supernovae is uncertain (e.g.. iBurrows et al. 
19951 : iMurphy et al.ll2004l ): if it is large enough, runaways produced by the supernova ejec- 
tion mechanism will also be single. For individual runaways, it is possible to identify a 
point of origin by measuring accurate distances, proper motions, and radial velocities (e.g., 
Hoogerwerf et al.l boOll : Ide Wit et all 1200,4 iMartinI bood: iHeber et al.ll2008h. For example, 
using Hipparcos data and accurate radial velocities, iHoogerwerf et al.l (120011 ) (i) associate 



the runaway star C Oph and the pulsar PSR J1932+1059 with a supernova in the Sco 0B2 
association and (ii) confirm that AE Aur, /i Col, and t Ori were ejected from the Trapezium 
cluster in the Orion nebula. Thus, both runaway mechanisms occur in nature. 

Existing observational and theoretical results place few constraints on the spatial and 
velocity distribution of runaways throughout the Milky Way. Because accurate proper 



in the solar neighborhood (e.g. 


Conlon et al. 


1990 




Holmgren et al. 


1992 




Mitchell et al. 


1998 


RoUeston e 


t al. 


1999; 


Hooeerwerf et al. 


2001 


: RamsDeck et al. 


2001; 


Maeee et al. 


2001; 


Lvnn et al. 


2004; 


Martin 


2004. 


2006 


) . Although some runaways have 


large distances from the 
, few are in the Galactic 


Galactic plane (e.g.. 


Greenstein & Sareent 


197^ 


; I 


ieber et al. 


2008) 



halo. Numerical simulati ons currently provid e little insight into runaways as probes of the 
Galactic halo structure. iDavies et al.l (120021 ) simulated the distri bution o f runa ways, but 
only in the context of high velocity white dwarfs near the Sun. iMartinI (120061 ) proposed 
using runaways to constrain the Galactic potential, but he concluded that known runaways 
provide few constraints on the large-scale properties of the Galaxy. 
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1.2. Prom the Galactic Disk to the Halo 



Recent observations of 'hypervelocity stars' (HVSs) in the halo motivate a broader 
investigation of the spatial and velocity distribution of runaways. The first HVS is a 3 
Mq main sequence st ar traveling at least twice the escape velocity of the Galaxy at its 



distance of ~ HQ kpc (IBrown et al. 



2005 . l2009al ). After two other serendipitous discoveries 



( iHirsch et al.ll2005l : lEdelmann et al.ll2005l ). subsequent targeted searches of the halo yielded 
a sample of 15 unbound HVSs with Galactic rest-frame velocities of 350-700 km s ~^ and a 



similar number of bound HVSs with rest-frame velocities of 275-350 km s ^ ( Brown et al. 



2006al lbl. [2007al lbl.[ 2009al ). Although the known HVSs are 3-4 Mq main sequence stars roughly 
uniformly distribute d in Galactic latitude, they are not isotropically distributed on the sky 
( IBrown et al. I l2009bh . 



The larg e spa c e velo city of HVSl suggests an origin in the Galactic Center. As first 
predicted by iHilld ( 1l988h , the tidal field of the massive black hole at the Galactic Center 
can unbind a close binary and eject one of the stars at velocities exceeding 2000 km s~^. To 
distinguish these high velocity stars from traditional runaway stars, Hills coined the term 
'hypervelocity star'. For typical ejection velocities expected from this mechanism (IHills 
Il988 ). we expect a range of velocities similar to those observed in H VSs. Although other dy- 
nami c al mechanisms involving a black hole can pro duce HVSs (e.g.. iHansen fc Milosavljevic 
20031 : IYu fc Tremaind l2003l : lO'Leary fc Loebl |2008| ) , the Hills mechanism makes clear pre- 
dictions for the exp ected velocity distribution of HVSs ejected from the Galactic Center. 
Kenyon et al.l (120081 ) use these predictions in a numerical simulation of the trajectories of 
HVSs through the Galaxy. They show that the relative number of bound and unbound HVSs 
predicted by the Hills mechanism agrees with observations of known HVSs. 

Recent observations suggest that runaways can also achieve unbound velocities. The 
likely unbound st ar HP 271791 wa s probably ejected from t he disk at 12-16 kpc from the 
Galactic Center (IHeber et al.ll2008l : iPrzybilla et al.ll2008al ). lJustham et al.l (120081) propose 



that the hot subdwarf US 708, with a heliocentric radial velocity of 708 km s ^ (IHirsch et al. 



20051 ). is also a runaway. The apparent overlap in the velocity and spatial distributions of 
runaways and HVSs suggests that several mechanisms may inject massive main sequence 
stars into the halo. The relative contributions of these processes to structure in the halo 
remains unknown. 

The distribution of runaw ays and HVSs is also an important is sue for large radial velocity 
surve ys. Surveys like RAVE JZwitter et al.ll2008h and SEGUE jAdelman-McCarthv et Z 



20081 ) measure radial velocities for hundreds of thousands of stars in the thin disk, thick 
disk, and halo. These data will provide fundamental constraints on the escape velocit y 
and total mass of the Milky Way JSmith et al.l I2OO7I : Ixue et al.l bood : Isiebert et al1l2008[ l. 
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Alth ough there are theo retical predictions for observable properties of HVSs in the halo 



[e.g., iKenyon et al.ll2008l ). there are no predictions for runaway stars in the halo. Thus, the 
contribution of runaways to velocity outliers in these surveys is unknown. 

Here, we use numerical simulations to make a first assay of the spatial and kinematic 
signature of runaways in large radial velocity surveys. Our focus is on intermediate mass 
1.5-6 Mq main sequence stars ejected from the disk into the halo. Intermediate mass stellar 
lifetimes are 10® - 10^ yr; thus, these stars formed recently in the disk. The supernova 
binary disruption scenario requires that the former companion of the runaway was a more 
massive star with a shorter lifetime. We do not consider massive (> 6 Mq) runaways, 
because they do not live long enough to reach large distances in the halo. We also do 
not consid er low mass stars because they are intrinsica lly faint and unobservable at large 
distances faoUmeier fc Gould! [20071 : iKenvon et ahlboosh . 



We construct a model for the spatial and velocity distributions of runaways. We use 
predicted ejection velocities from the supernova mechanism to eject stars from the exponen- 
tial disk and then we track their orbits through throughout the Galaxy. In §2, we describe 
the model and the results of our simulations. We apply the simulations to the runaway 
HD 271791 in §3, to HVS surveys in §4, and to halo radial velocity surveys in §5. We 
conclude in 56. 



2. The Simulations 



The model we construct to explore the global velocity and spatial distribution of run- 
aways in the Milky Way has two components: the gravitational potential of the Galaxy and 
the velocity distr ibution of ejected sta rs. For the Galactic potential, we use the 3-component 
model defined in lKenyon et al.l (120081 ). This disk, bulge, and halo model fits observations of 
the Galaxy on scal es from 5 pc to 1 00 kpc, and has a circular orbital velocity of 220 km s~^ 



at r = 8 kpc (e.g., iHogg et al.ll2005l ). 



We consider runaways ejected from an exponential disk with a representative distribu- 
tion of ejection velocities, explore the propagation of runaways in the Galaxy using a suite 
of simulations, and examine the resulting distributions in Galactocentric radius r and radial 
velocity Vrad at the end of the simulation. 
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2.1. Construction of the Simulations 



To generate populations of runaways in the Galaxy, we perform Monte Carlo simulations 
of 10^-10^ stars ejected into 3D orbits from the exponential disk. Each star begins on a 
circular orbit with velocity Wol'^mit) at some distance rinit from the Galactic center in the 
plane of the Galaxy. The star is ejected at an angle 9i relative to the orbital velocity vector 
and an angle hi relative to the plane of the Galactic disk. For a randomly oriented ejection 
velocity -Uej, the initial velocity of the ejected star is vinit = vq + Vej. For a star ejected in 
the direction of Galactic rotation {9i ^ 0), rjnit becomes the pericenter of its orbit; for a star 
ejected opposite Galactic rotation {9i ^ tt), r^nit becomes the apocenter of its orbit. 

To follow the trajectories of t hese ejected stars thr ough the Galaxy, we inte grate the 
equations of motion numerically. In lBromley et al.l (120061 ) and lKenyon et al.l (120081 ) . we used 
a simple leap-frog integrator to track stellar orbits through a ID Galactic potential. To pro- 
vide more accurate solutions for the trajectories of runaways in a 3D Galactic potential, here 
we use an adaptive fourth-order integrator with Richardson extrapolation (e.g., equations 
1-3 of iBromley fc Kenyon II2006I . see also Chapter 15 of Press et al. 1992). Starting at linit 
with velocity Vmit, the code integrates the full 3D orbit through the Galactic potential to 
track position and velocity as a function of time. We integrate the orbit for a random time 
t, where t lies between and tms, the main-sequence lifetime of the runaway star. For stars 
ejected during the supernova explosion of a more massive and much shorter-lived companion 
star, t is roughly the age of the runaway star. 

To motivate choices for the initial conditions of our simulations, we derive the approx- 
imate velocities of runaways capable of reaching the outer halo of the Galaxy. Numerical 
simulations of binary systems disrupted by su pernova explosions show that runaways have 
maximum ejection velocities of ~ 400 km s~^ (jPortegies Zwartll2000l : iDrav et al.l 12005). Re- 



alistic dynamical ej ection models yield similar maximum ejection velocities (ILeonard 



1993 



Fregeau et al.ll2004j ). If the highest velocity runaways are ejected in the direction of Galactic 
rotation, they receive an additional kick of ~ 200 km s~^. Thus, the maximum possible 
ejection velocity is roughly Vej^max ~ 600 km s~^. Unb ound HVSs in the outer halo have 
Vrad ^ 400 km s~^ at r ^ 70 kpc (IBrown et al.ll2009al ). To match these properties in our 
model for the gravitational potential of the Milky Way, runaways must have initial velocities 
Vej = 725 km s 



at rinit= 1 kpc, Vpj = 550 km s ^ at rinit= 10 kpc, and v, 



500 km s"^ at 

30 kpc (see Fig. 2 of iKenyon et al.ll2008l ). Thus, runaways ejected from r^nit < 5 kpc 



cannot reach the outer halo. This conclusion leads us to consider a fiducial set of simulations 
with Vcj= 400 km s~^ and rjnit = 10 kpc. 

Long travel times through the halo constrain our choices for the initial ages and masses 
of main sequence runaway stars. Stars with Vej = 400 km s~^ travel ~ 40 kpc in ~ 100 Myr. 
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Thus, stars with main se quence lifetimes tm^ > 100 Myr and ini t ial masses m < 4-5 M (7, 
can reach the outer halo (ISchaller et al.lll992l : ISchaerer et al.lll993l : iDemarque et al. 1 120041 ) . 
Stars with m ~ 4-6 Mq and tm „ ^ 65-160 Myr can reach the inner halo. To avoid confusion 
with indigenous halo stars (see lBrown et al.ll2006al . l2009al). we focus on main sequence stars 
w ith m = 1.5-6 an d tms = 65 Myr to 2.9 Gyr JSchaller et al.lll992l : ISchaerer et aDll993 



Demarque et al. II2004I . see also Table 1). These lifetimes are much longer than the lifetimes, 



t 



m.9 5, 10-20 Myr, o f their m > 10 MfT^ binary companions that explode as supernovae 
( ISchaller et al.lll992l : ISchaerer et al.l 119931 ). Thus, we assume for simplicity that runaw ays 
have ages of zero when they are ejected from the binary (see also iPortegies Zwartll2000l ). 



2.2. Results for a Fiducial Simulation 

We first investigate the Galactic distribution of runaways for the fiducial parameters 
rinit = 10 kpc, Vej = 400 km s~^ and m = 3 Mq. Figure [T] shows results for 10^ runaways; we 
plot Vrad as a function of r. In the upper panel, we color-code the points in four intervals of 
9i for stars with bi < 30°. Blue (yellow) points indicate stars ejected with 9i = — 7r/4 ^ 7r/4 
{6i = Stt/A —>■ 57r/4). Cyan and magenta points indicate stars ejected at intermediate angles, 
6i = ±(71/4 7i'/2) for cyan points and 6i = ±(vr/2 37r/4) for magenta points. In the 
lower panel, we color-code the points in three intervals of bf bi < 30° (cyan), 30° < bi < 
50° (magenta), and bi > 50° (yellow). 

The distribution of runaway stars in Fig. [1] has several features. In the top panel, stars 
with the largest initial velocities (blue points) are marginally bound to the Galaxy. Thus, 
they travel far into the halo, reach a maximum distance r^ax, and then fall back towards the 
Galactic Center. However, the travel time for this orbit, t ^ it r^ax / Vmit ~ 1 Gyr, is longer 
than the main sequence lifetime of a 3 M© star, tms ~ 350 Myr. These stars can only be 
observed on their outward path through the halo and always have positive radial velocities. 
Thus, the sequence of blue points from (r, v^ad) = (20, 500) to (160, 340) is an age sequence, 
with 20-30 Myr-old stars at small r and 300-350 Myr-old stars at the largest r. 

Stars ejected opposite to the direction of Galactic rotation (yellow points) have the 
smallest initial velocities and thus cannot reach large r. These stars live long enough to 
orbit the Galaxy at least once. Thus, the group of yellow points represents a mixture of 
young and old stars with a velocity distribution symmetric about zero. 

Stars ejected at intermediate angles (cyan and magenta points) fill the (r, Vrad) space in 
between the group ejected along or against Galactic rotation. For each group, the sequence 
from largest Vrad to largest r is an age sequence, with younger stars at small r and older stars 
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at large r. Thus, most of the lower envelope of the complete ensemble of points, extending 
from (r, Vra.d) = (30, —200) to (r, Vj-ad) = (160, 340) consists of old stars near the end of 
their main sequence lifetimes. 

The lower panel of Fig. [1] illustrates how the maximum radial velocities decrease with 
increasing bi. Stars at high latitude (yellow points; bi > 50°) have the smallest velocities, 
with Vrad ^ 400 km s~^ at all r. Stars at intermediate latitudes (magenta points; 30° < bi < 
50°) reach Vrad ^ 450 km s~^ at all r; stars at low latitudes (cyan points; bi < 30°) have 
the largest final velocities (vrad ~ 500 km s~^) and reach the largest final radii, r = 150- 
160 kpc in 350 Myr. The form of the Galactic potential produces this variation. Stars ejected 
perpendicular to the Galactic disk receive a smaller kick from Galactic rotation compared 
to stars ejected parallel to the disk. Once stars leave the plane of the disk, they also feel the 
full disk potential. Thus, stars ejected perpendicular to the disk have smaller initial radial 
velocities and decelerate faster than stars ejected parallel to the disk. 

To show how the results depend on each input parameter. Figure |2] plots the velocity 
distributions for simulations where we vary each of the three fiducial parameters separately 
and hold the other two fixed. Here, we color code the points with their final galactic latitude. 
Instead of illustrating relative densities, our intent in this figure is to show changes in the 
shape of the distributions of Vrad and r as functions of the various input parameters and bj, 
the Galactic latitude of stars at the end of the simulation. The upper left panel of Figure [2] 
repeats the distribution of points in the lower panel of Fig. [T] with a different color coding. 
Here, the variation of v^ad and r^ax with final latitude is much more pronounced than the 
variation with ejection angle. 

The predicted distributions clearly depend on rinit and Vinit (upper right and lower left 
panels of Figure [2]). Stars ejected at small r have to climb out of a deeper potential well than 
stars ejected at large r. At fixed Vej, runaways ejected from 30 kpc have larger velocities 
in the outer halo (maximum Vrad = 500 km s~^) than stars ejected from 10 kpc (maximum 
Vrad = 400 km s~^). Similarly, stars with smaller (larger) ejection velocities maintain smaller 
(larger) velocities throughout their passage through the Galaxy. Thus, at fixed rinit, runaways 
ejected at larger velocities have larger velocities in the outer halo and reach larger distances 
in the halo. 

Finally, stellar lifetimes have a significant impact on the distributions of position and 
velocity (Figure [21 lower right panel). For fixed initial ejection velocities, stars with longer 
main sequence lifetimes travel farther into the halo. Stars with longer lifetimes can reach 
large distances with smaller initial radial velocities. Thus, lower mass stars have smaller v^ad 
at larger r. For our adopted main sequence lifetimes, 1.5-4 runaway stars ejected from 
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10 kpc have asymptotic radial velocities 

t^rad(w,io) ~ 250 {MJ2 M^f'^ km s"^ (1) 

at maximum Galactocentric distances 

Tmax^io ~ 400 (M,/2 M0)-2-4 kpc . (2) 

At 30 kpc, runaways have larger asymptotic radial velocities 

v,^{r^a.m) - 400 {MJ2 M^f'^ km s"^ (3) 

at maximum Galactocentric distances 

~ 575 (M,/2 M^Y^-^ kpc . (4) 

In both equations for Tmax-, the l arge exponent follows from the mass dependence of the main 
seque nce lifetime {tms oc M~^] ISchaller et al.lll992l ; ISchaerer et al.l Il993l ; iDemarque et al 



20M). 



Tables [TH2] quantify these generic conclusions for an expanded set of fiducial simulations. 
For simulations with rinit = 10 kpc, the median velocities in Table [1] show clear trends with 
stellar mass and r. Only the highest velocity runaways can reach large r; for all stellar masses, 
the median velocity increases with increasing r. Because their main sequence lifetimes are 
longer, lower mass stars to orbit the Galactic Center 2-3 times. Thus, median velocities for 
lower mass stars are closer to zero than those for higher mass stars. 

The trends of median radial velocity with hf depend on r and stellar mass. Stars 
ejected into the disk have larger initial velocities and always travel farther than stars ejected 
perpendicular to the disk. Thus, at large r, stars at smaller galactic latitude have larger 
median Vj-ad- Higher mass stars show the largest variation of median Vj-ad- For every 30° 
increase in 6/, the median Vrad declines ~ 25-75 km s~^ for 1.5-3 stars. At small r, 
1.5-3 stars with hf < 30° and bf > 50° were ejected against Galactic rotation; these 
stars have a range of orbital phases and thus have median v^ad close to zero. At intermediate 
bf, there is a mixture of stars ejected opposite to Galactic rotation (which have median Vrad 
close to zero) and stars ejected with Galactic rotation (which have larger median Vrad)- Thus, 
the median velocities for bf = 30°-50° are generally large. 

The typical radial velocity dispersions of runaway stars also vary consistently with r, 
bf, and stellar mass (Table [2]). Stars observed at small r are a mix of young stars ejected 
into the outer Galaxy and older stars orbiting in the inner Galaxy. Thus, these stars have 
large velocity dispersions for all stellar masses. Because massive stars have shorter lifetimes. 
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they do not travel far from their ejection point and tend to have smaller velocity dispersions. 
Stars observed at large r are older stars traveling on extended bound orbits. None of these 
stars live long enough to fall back into the Galaxy. Thus, all have large, positive Vj-ad and 
small velocity dispersions. 



2.3. Results for a Galactic Distribution of Runaways 



We now consider a simulation for a complete ensemble of runaways ejected from the full 
Galactic disk. To make this simulation, we adopt probability functions to assign rinit and Vgj 
for each runaway star, integrate the orbit through the Galaxy, and derive the radial velocity 
Vrad and position r for e ach runaway at a random time t in its orbit. Our approach differs 
from lDavies et al.l (120021 ). who simulated runaways ejected uniformly from 1 < r < 10 kpc. 



For the initial distribution of rj^j t,, we assume stars are ejected from an exponential disk 
with a radial scale length of 2.4 kpc (Siegel et al. 20021). We adopt an inner radius of 3 kpc 
and an outer radius of 30 kpc ( Brand &: WouterlootI 2007 ). Thus, the probability distribution 
for rinit is 

p(ri„it) oc r,„,i e-^-*/2.4 kpc 

Runaways ejected from rjnit < 3 kpc cannot reach the outer halo. Runaways ejected from 
Tinit > 30 kpc are too rare to make a significant impact on the derived distribution of (r, 

Vrad)- 

For the initial distribution of Vej , we adopt results from published analyses of runaways 
ejected from binary systems disrupted by supernovae. Although dynamical encounters can 
also produce runaways, there are no published simulations predicting p(wej) the probability 
distribution of ejection velocities. For Vej ~ 20-400 km s~^, a simple function, 



p{Vc 



OC e 



-Ucj/150 km s 



(6) 



provides a reasonable match to the iPortegies ZwartI (120001 ) simulations of binary supernova 
ejections. 



We select the IPortegies ZwartI (120001 ) ejection velocity distribution function because 
it is physically well motivated and because a similar distribution function is not available 



for the dynamical e 



( iLeonard fc Duncan 



ection mechanism. Dynarnical ejections can attain higher velocities 



19881 . Il990l : iLeonardI Il99ll . Il993l ): however, the theoretical maximum 



ejection velocity is not realizable b ecause compact bin ary interactions are more likely to 
merge stars than to eject runaways (jFregeau et al.ll2004j ). Moreover, the ejection rate from 
binary- binary encounters is probably smaller than supernova ejections for intermediate mass 
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stars. Dynamical ejections depend o n the joint probabili ty of colliding two binaries within the 



main sequence lifetime of the stars (IBrown et al.ll2009al ). Thus we use the binary-supernova 



ejection velocity distribution as representative of the runaway process for 1.5-6 Mq stars. 

Full-disk simulations of runaways yield many of the same features in the Vrad-^ diagram 
(Figure [3], top panels). As in Figure [21 the lower envelope of the set of points is defined by 
an ensemble of old runaways close to the end of their main sequence lifetimes. Stars ejected 
into the disk (cyan points) receive the maximum ejection velocity. Thus, these stars have 
the largest Vrad at all r. Stars ejected into the halo (yellow points) have the smallest ejection 
velocity, the smallest v^ad at all r, and the smallest radial extent. More massive stars with 
shorter tms do not live long enough to reach large r. Thus, lower mass stars have more 
extended radial distributions than more massive stars. 

The major difference between the fiducial and the full-disk simulations is the variation 
of the maximum Vrad with r. Runaways ejected from the inner disk slowly decelerate as they 
move from their point of origin. In an ensemble of runaways ejected from a single radius in 
the inner disk, stars at larger r therefore have smaller v^ad than stars at smaller r (upper left 
panel of Figure [2]). However, runaways ejected from the outer disk coast outward at roughly 
constant Vrad (upper right panel of Figure [2]). These two features of the evolution combine 
to produce a roughly constant maximum Vrad with r in the top panels of Figure [3l This 
maximum velocity is independent of stellar mass and is roughly 100 km s~^ larger than the 
maximum ejection velocity. 

Tables |3H5] list the median velocities, radial velocity dispersions, and vertical velocity 
dispersions for the full-disk simulations of 1.5-6 Mq stars. For all stellar masses, stars ejected 
at larger Vrad reach larger r. Thus, the median radial velocity increases with r. Because the 
most distant stars must have roughly the same high ejection velocity to reach large r, these 
stars have smaller velocity dispersions than the mix of low and high velocity runaways at 
small r. For stars ejected with similar velocities, lower mass stars live longer and can reach 
larger r. Thus, the median v^ad and the velocity dispersions increase with increasing M*. 

The distribution of runaways in the Vrad-''^ diagram differs from the distribution of HVSs 
(Figure [31 lower right panel). The lower envelope of the HVS distribution is identical to the 
runaway distribution and is composed of stars with ages comparable to their main sequence 
lifetimes. However, the distribution of HVSs has three features not observed in runaway stars. 
The HVS distribution has a core of stars at r < 3 kpc; these stars were eje cted from the 



Gala ctic Center at small velocities (< 700 km s ^) and cannot reach large r ( iKenyon et al. 



20081 ). Because some HVSs are ejected at very high speeds (> 1200 km s~^), these stars 
can reach large r with radial velocities much larger than any runaway star. Finally, the 
velocities of HVSs ejected isotropically from the Galactic Center are independent of b. Thus, 
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high velocity HVSs are observable at all b. In contrast, high velocity runaway stars are only 
observable at low galactic latitude. 

The velocity distribution of runaways is also very different from the velocities of halo 
stars (F igure [Si lower left panel). Halo st ars have a ID velocity dispersion of rou ghly 100-110 
km s ~^ (jHelmill2008l : iBrown et al.ll2009al ) and a spatial density close to n oc (jSiegel et al. 
20021 ). Thus, the envelope of the halo velocity distribution is symmetric about zero and 
declines slowly with radius. At large r, there are many halo stars with v^ad ~ 0. In contrast, 
there are no intermediate mass HVSs or runaways with v^ad ~ at large r. 

To quantify the differences between runaways, HVSs, and halo stars, we now consider 
the predicted radial surface density of runaways (Figure H]). Compared to the initial surface 
density of the Galactic disk (dot-dashed line in Figure H]), runaways are much more radially 
extended. Runaways ejected against Galactic rotation populate the inner disk (r < 3 kpc). 
These runaways have a large velocity dispersion at all b (top panels of Figure [3]). Runaways 
ejected at high velocity along Galactic rotation populate the outer disk. These ejections 
produce a power-law surface density profile, E oc r^" with n = 3.5-3.6, at intermediate r 
and an exponential decline at large r. Longer-lived lower mass stars have the most extended 
power-law component. Short-lived massive stars have exponential density profiles more 
similar to the density profile of the Galactic disk. 

Although the density profiles of 1.5-3 Mq runaways are extended, they are much steeper 
than the density profiles of halo stars and HVSs. Fits to obser vations of halo stars typically 
yield power law density profiles with n ^ 2.7-3.5 (jHelmill2008l . and references therein). Our 
numerical simulations suggest that low mass HVSs have a boun d component with n = 3 
and an unbound component with n = 2-2.5 (IKenyon et al.ll2008l . see also Hills 1988). The 
short, finite lifetimes of massive stars steepens the HVS density profile ('^ ^ 3) at r > 80 
kpc. However, these density profiles remain shallower than the density profiles of runaway 
stars. 



The predicted vertical density distribution also distinguishes runaway stars from halo 
stars and HVSs (Figure E]). This density distribution has two components. Low velocity 
runaways produce a 'thick disk' with a vertical scale height of 300-1000 pc. Although more 
distant runaways have slightly larger 'thick disk' scale heights, the scale height is independent 
of stellar mass. High velocity runaways lie in an extended disk-shaped halo with a vertical 
scale height of 2-40 kpc. More distant runaways have much larger vertical scale heights. At 
r < 50-60 kpc, the scale height of the extended halo is independent of stellar mass (Table E]). 
At r > 60 kpc, the vertical scale height depends on stellar mass. Massive stars (M^ > 4 Mq) 
will short main sequence lifetimes cannot reach large r and thus have no measurable scale 
height. For stars with longer stellar lifetimes, lower mass stars have smaller scale heights at 
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large r. These smaller scale heights result from low velocity ejected stars, which have time 
to reach large r only for the lowest mass main sequence stars. 

The disk-shaped density distribution of the highest velocity runaways differs from the 
spherically symmetric density distributions of halo stars and HVSs. In our simulations, 
HVSs ejected from the Gal actic Center have a spherically symmetric, power law density 
distribution with n ^ 2-2.5 (iKenyon et al.ll2008l ). Halo stars are also distributed s pherically 
symmetrically and have a steeper radial density profile {n ^ 2.7-3.5: iHelmil 120081 ). Because 
HVSs have a shallower density profile than halo st ars, it is easier to identify HVSs at large 



halo distances than at small halo distances (e.g., iBrown et al.l l2005l : iBromley et al.l 12006 



Kenyon et al.ll2008l ). The steeper density profiles produced in our runaway star simulations 
suggest nearby runaways are easier to identify than distant runaways. We consider this 
possibility further in §4-5. 

To conclude this section, we examine several additional properties of simulated runaways 
in the outer Galaxy. Most runaways with r > 60 kpc are ejected from inside the solar circle 
(Figure [H]). Although most ejected stars have rjnit ~ 3-6 kpc, small ejection velocities prevent 
them from reaching r > 60 kpc (§2.1). Because the depth of the Galactic potential is smaller 
for runaways with rinit ~ 10-20 kpc, these stars make up a large percentage of stars with 
r > 60 kpc (see also Figure [2]). In addition, most runaways at large r are low mass stars 
with stellar lifetimes long enough to reach the outer Galaxy. With median radial velocities 
of 10-200 km s~^ (Table [3]), nearly all of these runaways are bound to the Galaxy. Before 
reaching the outer galaxy, massive stars {M^ > 4 Mq) evolve off the main sequence and are 
unobservable at r > 60 kpc. 

Our simulations yield a small fraction - 0.07% - of runaways that are not bound to the 
Galaxy. More than half of unbound runaways are ejected from outside the solar circle (double 
dot-dashed line in Figure E]). Nearly all unbound runaways have low Galactic latitude (6/ < 
30°; Figure [7] and Table [7]). For 1.5-6 Mq stars, the shape of the cumulative probability 
function for 6/ is nearly independent of stellar mass. The median of the distribution, however, 
varies slowly with stellar mass. Our results suggest median Galactic latitude bf^med = 9°-10° 
for 1.5-3 Mq stars and bf med = 7° for 6 Mq stars (Table [7j). 

These results for the distribution of bj_ contrast with numerical simulations of HVSs, 



which yield a uniform distribution in bf ( iBromley et al.l l2006l : iKenyon et al.l |2008| ) . In the 
Hillsl ( 1l988l ) ejection mechanism, HVSs are ejected isotropically from the Galactic Center. 
For these stars, the disk is a small perturbation on the potential. Thus, simulated HVSs 
remain uniformly distributed in 6j throughout their path through the Galaxy. For runaway 
stars. Galactic rotation provides a significant kick to the ejection velocity. Stars ejected with 
this rotation are more likely to have unbound velocities than other stars (Figure [TH3]). Thus, 
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unbound runaways are confined to the plane of tlie disk. 

Altliougfi unbound runaways liave low bj independent of M^, the predicted distances 
are sensitive to (Figure [Hand Table [7]). The results in Figure [S] show that the cumulative 
probability of r grows roughly linearly from r ^ 6 kpc to r 0.75rmax and then asymp- 
totically approaches 1 at r ^ Tmax- The limiting distance Vmax of unbound runaways is a 
strong function of stellar mass. Long-lived low mass stars travel far into the outer Galaxy 
{"^max ~ 500-1000 kpc for M^= 1.5-2 Mq); short-lived high mass stars rapidly evolve off the 
main-sequence and cannot reach large r {vmax ~ 50-70 kpc for = 5-6 Mq). 

2.4. Summary of the Simulations 

Our simulations demonstrate how stellar evolution and the Galactic potential combine to 
influence the dynamical properties of runaway stars ejected from the Galactic disk. Runaways 
that receive the maximum kick from the binary-supernova mechanism, ^ 400 km s~^, can 
travel from the disk into the halo. These stars produce an extended disk-shaped distribution 
of stars, where the radial and vertical scale lengths are much larger than those of the main 
stellar disk. The size of this extended disk is very sensitive to stellar mass. Massive stars 
with short stellar lifetimes are much less extended than long-lived low mass stars. Because 
runaways ejected along the direction of Galactic rotation have higher ejection velocities and 
climb out of a shallower potential well than other runaways, these stars reach larger distances 
in the outer Galaxy. Although high velocity runaways appear at all bf, the fastest unbound 
runaways are at low galactic latitude, bf < 30°. 

Comparisons of our results with observations of halo stars and simulations of HVSs 
suggest clear differences between the three populations. Halo stars and HVSs are uniformly 
distributed in b; runaways are concentrated in the disk. The radial density gradients of halo 
stars and HVSs are shallower than those of runaways. The Vj-ad distributions of HVSs and 
runaways are clearly non-gaussian compared to observations of halo stars; HVSs tend to 
have larger Vj-ad than runaways. These differences suggest clear observational discriminants 
of the populations, which we explore in the following sections. 

3. Application to the Hyper-Runaway HD 271791 



For a first application of our simulations, we consider the runaway B star HD 271791. 
This B2-3 HI star lies well below the Galactic plane jz ~ —10 kpc) and has a large he - 
liocentric radial velocity of 442 km s~^ (IKilkenny fc Stondll988l : [Kilkenny fc Mullerlll989l ). 
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The observed proper motions suggest a large range in Galac tic rest-franie velo city, 530- 
920 km s-\ and an origin in the outer disk at r ^ 12-16 kpc (iHeber et al.ll2008l ). Detailed 
abundance analyses yield a subsolar iron abundan ce and enhanced abundances of the a- 
capture nuclei O, Ne, and S (jPrzybilla et al.ll2008al ). The subsolar iron abundance supports 
an origin in the outer Galaxy; the high abundances of a-nuclei imply contamination of the 
atmosphere from a nearby supernova. 



There are two proposals for the origin of HD 271791 as a hyper-runaway star. iPrzybilla et al 



(l2008al ) suggest that HD 271791 is the secondary of a massive binary disrupted by a super- 
nova. However, the upper limit of the obser ved rest-frana e veloc ity, ~ 900 km s~^, is hard 
to achieve in a close binary system. Thus iGvaramadzd (120091 ) prefers ejection during a 
three-body or four-body encounter in the dense core of a massive star cluster. 

The large range of allowed rest-frame velocities does not permit a unique interpretation 
for this star. Here, we first consider the minimum rest-frame velocity of 530 km s~^ and note 
how larger velocities impact our conclusions. 

For a rest-frame velocity of 530 km s~\ HD 271791 is a plausible runaway star produced 
by a supernova explosion in a massive binary system. With a distance r 21 kpc from th e 
Galactic Center and a height z ^ —10 kpc below t he Galactic pla i ie (jHeber et al.ll2008l ). 
the star is marginally bound to the Galaxy using the iKenyon et al. torn potential model 
for the Milky Way. The radial velocity is similar to the maximum velocity achieved in our 
simulations (Figure [3]). The galactic latitude of 6 ~ —25° is marginally consistent with the 
runaway interpretation. Because the highest velocity runaways are ejected into the disk, 
only ^ 5%-10% of unbound runaways have bf > 20° (Figure [7]). Roughly 30% of unbound 
runaways are ejected from the apparent origin of HD 271791 at r ^ 12-16 kpc (Figure [6]). 
After the ejection, an 11 Mf^ star can reach z 10 kpc during its main sequence lifetime 
of 20-30 Myr (Table 3, and iHeber et al.ll2008l ). Thus, the observed spectral type is also 
consistent with the runaway interpretation. 

Larger rest-frame velocities weaken the case for the supernova ejection model. In our 
simulations, the largest observed rest-frame radial velocity is roughly 100 km s~^ larger 
than the maximum ejection velocity from the binary. An observed rest-frame velocity of ~ 
900 km s~^ requires an ejection velocity of 800 k m s~^, roughly a facto r of two larger than 
the maximum achieved in numerical simulations (jPortegies Zwartll2000l . see also Przybilla et 
al 2008). Off-center supernova explosions might enhance the ej ection velocity, bu t factor of 
two increases for a massive star like HD 271791 seem unlikely (IGvaramadzd 120091 ) . 



Dynami c al ejec tion mechanisms can explain the high rest-frame velocity of HD 271791. 
Gvaramadzd (120091 ) outlines several mechanisms where the interactions among binary or 
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triple systems produce an ejected star with a velocity of < 800 km s~^. To account for the 
enhanced a-nuclei in the atmosphere, the dynamical encounter(s) in each mechanism must 
occur close in time to the supernova explosion of one of the companion stars to HD 271791. 



Low prob abilities complic a te all formation mechanisms for HD 271791 (IBrown et al. 



2009al ). In the iPrzybilla et al.l (j2008al ) model, a supernova explosion in a 80 M© star pro- 
duces the high ejection velocity of the secondary. For a Galactic star formation rate of 4 
Mq yr~^ (e.g., IPiehl et al.ll2006l ). we expect ^ 5 x 10"^ stars with > 80 Mq during the 
25 Myr ma in sequence life time of HD 271791. Roughly one-third of all O star binaries are 
twins; for alSalpeterl ( 1955 ) initial mass function, 5% of the rest have 10-12 M0 companions 
(iKobulnicky fc Fryerl 120071 ). Thus, ~ 1000 massive binaries with 80 Mq primary stars have 
10-12 M0 secondaries. If all supernova ejections in these binaries produce a runaway, we 
expect 1000 x 0.08% ~ 1 unbound runaway similar to HD 271791 every 25 Myr. Extrapo- 
lating the results of our simulations to 11 Mq runaways, the joint probability of observing 
this runaway at 6 ^ —25° and at the end of its main sequence lifetime at r 20 kpc is ~ 
0.1 (see Table [7]). Thus, HD 271791 is an unlikely runaway star. 



Brown et al.l (j2009al ) show that the likelihood of observing a hyper-runaway from a 



dynamical ejection is also very small. Dynamical ejections require interactions between two 
binary systems composed of massive stars. Thus, the probability of a dynamical ejection is 
the joint probability of interacting masive pairs of binaries within the stars' main sequence 
lifetimes. In an ensemble of dense clusters capable of producing high velocity ejections, the 
probability of a 3-4 Mq hyper-runaway is ~ 10~^. For a Salpeter (1955) initial mass function, 
11 Mq stars are ~ 4-5 times less likely than 3-4 Mq stars. Thus, dynamical ejections are 
much less likely than supernova ejections. 

An improved proper motion for HD 271791 would place better constraints on formation 
mechanis ms. Current data have large uncertainties, leading to a large range in rest-frame 
velocity (IHeber et al.ll2008l ). Observations with GAIA, scheduled for a 2011 launch, would 
yield a very accurate rest-frame velocity. 



4. Application to Hypervelocity Stars 

We now explore the observational consequences of runaway star distributions for HVS 
surveys. First, we compare our runaway simulations to numerical simulations of the apparent 
magnitude and heliocentric radial velocity distributions of HVSs. We then compare our 
runaway simulations with the HVS observations of Brown et al. Finally, we estimate the 
possible contribution of runaways to the observed sample of HVSs. 
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4.1. Predicting Observables from Simulations 



We begin by "observing" our numerical simulations for 3 runaways and HVSs in 
a heliocentric reference fr ame. K n own HV Ss come mostly from the radial velocity survey 
of brown et all J2006al jbl. l2007al lbl. l2009ah who target objects with th e colors of 3-4 
stars. Three of the HVSs are confirmed ~3 M f:^ main se quence stars (IFuentes et al.l 12006 



Przvbilla et aPboOSbl : llopez-Morales fc BonanSboOsh . Icirardi et all J2002I . booj ) stellar 
evolutionary tracks show that a solar metallicity, 3 M0 star spends 350 Myr on the main 
sequence with an average luminosity of Mg = 0.0. 

Thus we calculate heliocentric distances and apparent magnitudes for the simulated 
runaways assuming Mg = 0.0 appropriate for a 3 M0 star. We assume the Sun is located at 
r = 8 kpc. We shift the origin of the coordinate system to the Sun, and derive heliocentric 
radial velocities by taking v ■ r. 

Figure plots a heliocentric view of the simulations as a function of Galactic longitude, 
latitude, and apparent magnitude. The top panels of Figure [9] plot the number distribution 
of runaways and HVSs. The bottom panels of Figure [9] plot the 50*^^, 90**^, and 99*^^ percentile 
heliocentric radial velocity of runaways and HVSs. 

The spatial concentration of runaways in the disk, discussed in §2 above, is evident in 
Figure m Our models predict a greater fraction of runaways in the Galactic center hemisphere 
|/| < ±90° and at low Galactic latitudes b < 30°. The fraction of HVSs, on the other hand, 
is larger than the fraction of runaways in the direction of the Galactic anti-center / = 180° 
and the Galactic pole |6| = 90°. Compared to HVSs, runaways are also apparently bright: 
80% of 3 Mq runaways are brighter than g = 16, whereas 85% of 3 Mq HVSs are fainter 
than g = 16. 

Runaways are ejected from a rotating disk. This rotation is apparent in the distribution 
of heliocentric radial velocities: the median (50*^ percentile) runaway velocity is negative 
in the direction of Galactic rotation 0° < / < 180° and positive in the direction opposite 
Galactic rotation 180° < / < 360°. Similarly, the latitude dependence of runaway velocities 
reflects the boost from Galactic rotation at low latitudes. HVSs, on the other hand, are 
ejected on purely radial trajectories and show no rotation. Median HVS velocities exceed 
the 99*^ percentile runaway velocity in every direction on the sky. Fainter stars are faster 
because only the fastest runaways and HVSs survive to reach the largest distances. 
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4.2. Comparing Simulations and Observations 



We now compare our runaway simulation of 3 M© stars to observations of HVSs. Ob- 
served HVSs are significant velocity outliers with mi nimum radial velocit ies in the Galactic 
rest frame > +400 km s~^. The well-defined survey of iBrown et al. fl2006al lbl. l2007al jbl. boOQah 
samples stars with magnitudes 15 < g < 20.5 and latitudes 30° ^ b < 90°. 

Over the range of magnitude and latitude sampled by the Brown et al survey, only 
0.001% of the simulated 3 Mq runaways have velocities > +400 km s~^. These runaways 
differ in three ways from the observed HVSs. (1) Simulated runaways with radial velocities 
exceeding 400 km s~^ are located at low latitude with median b = 34 °; observed HVSs ar e 
distributed uniformly across Galactic latitude with median b = 51° (IBrown et al.l l2009bl ). 
(2) The > 400 km s"^ simulated runaways are brigh t with median g = 16; observed HVSs 
are faint with median = 19 (IBrown et al.l l2009a] ). (3) The fastest simulated runaway 
with b > 30° has velocity +450 km s~^; observed HVSs have velocities up to +700 km s~^. 
We conclude that runaways cannot significantly contaminate the HVS samples because the 
observed distribution of HVSs differs so markedly from that expected for runaways (Figure 

We now consider how many runaways with radial velocities exceeding +400 km s~^ 
might be included in the Brown et al HVS survey. To explore this point, we normalize the 
number of runaways in our simulation to the total number of 3 M0 stars forme d in the last 



350 Myr. Assuming the star formation rate in the Galactic disk is 4 M0 yr ^ (jPiehl et al. 
20061), 1.4 X 10^ M0 of stars have formed in the disk in the past 350 Myr. A Salpeter initial 



mass function, integrated from 0.1 to 100 Mq and normalized to 1.4 x 10^ Mq, predicts 
1.3 X 10'^ 3-4 Mq stars. Assuming that ~1% of all stars are ejected as runaways (see 
below), we predict ~1 3 Mq runaway with +400 km s~^ radial veloci ty in the Brown et a l. 
HVS survey. This prediction is comparable to the analytic estimate in lBrown et al.l (j2009al ). 
and suggests that perhaps one of the observed HVSs may be a runaway. 



5. Application to Halo Radial Velocity Surveys 



Large modern photometric and radial velocity surveys open the possibility of global 
constraints on the number and distribution of runaways. In this section, we consider several 
examples. We first compare the predictions of our runaway star simulations with the stellar 
population of the Milky Way derived from color selected star counts in the SDSS. This com- 
parison indicates where future searches for ru naways would be m ost productive. Combined 
with a 2MASS-selected spectroscopic survey (IBrown et al.l 120081 ) . the SDSS counts provide 



- 19 - 



an upper limit on the fraction of disk A main sequence stars ejected as runaways. 

Our numerical simulations demonstrate that runaways reach the inner halo of the Milky 
Way. Their velocity dispersion is comparable to the halo velocity dispersion; their rotation 
is comparable with the rotation of the thick disk. Thus runaways are difficult to identify by 
their kinematics, but their high metallicities should distinguish them from typical halo stars. 



5.1. Comparing Runaways to Star Counts 



We begin by comparing our runaway simulations to the observed stellar population of 
the Milky Way as revealed by star counts in th e Sloan Digital Sky Survey (SDSS) Data Re- 
lease 6 (DR6) (lAdelman-McCarthy et al.ll2008l ). For this comparison, we select simulations 
for 1.5 Mq, 2 M0, and 3 M0 stars. As before, we "observe" the simulations from a helio- 
centric reference frame. We calculate apparent magnitudes assuming the 1.5 Mq, 2 Mq, and 
3 M^ runaways have main sequence luminosities of Mg = +2.9, +1.5, and 0.0, respectively 
dCirardi et allbood l2004[ ). 



To compare observed number counts of stars with predictions from our numerical sim- 
ulations, we consider only runaways that fall in the region of sky imaged by the SDSS 
DR6. For the stars in the SDSS, we count those stars with colors 0.15 < {g — r)o < 0.20, 
-0.15 < - r)o < -0.05, and -0.35 < {g - r)n < -0.25 appropriat e for 1.5 Mq, 2 Mq, 



and 3 Mq stars, respectively, according to the iGirardi et al.l (|2002| . |200J) stellar evolutionary 
tracks for solar metallicity stars. 

Figure [TU] plots the resulting number fraction of runaways and SDSS stars as a function 
of Galactic latitude and apparent magnitude. Over the region surveyed by SDSS, a larger 
fraction of runaways are found at low latitudes compared to star counts. This latitude 
dependence reflects the flattened distribution of runaways compared to the population of 
halo stars that dominate the star counts. Runaways contribute a negligible amount to the 
observed stellar (halo) population at faint magnitudes, g ^ 17. 



5.2. Upper Limit on the Runaway Fraction 



The runaway fraction of O- and B-type stars ha s long been known to be ~40% and 
■'5%, respectively (lBlaauwlll96ll : iGies fc BoltonI 119861). However, there are few comparable 



constraints on the fraction of A-type runaways. IStetsod (Il98ll ) used tangential velocity 
cuts of bright stars with rough spectral types to estimate that at most ~0.3% of solar 
neighborhood A-type stars are runaways. Here, we combine constraints from our simulations 
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and existing spectroscopic surveys to place an independnent upper limit on the fraction of 
2 Mq stars ejected as runaways. 



To make this estimate, we use the iBrown et al.l (12008) spec trosco pic survey of A-type 
stars in the Two Micron All Sky Survey (jSkrutskie et al.ll2006l ). The iBrown et al.l (120081 ) 
survey is complete over 4300 deg^ to a mag nitude limit of Jo = 15.5, equivalent to ~ 15.5 
for a zero-color A star. iBrown et al.l (120081) find that 40% of A-type stars at 15**^ mag are 
main sequence stars (the other 60% are evolved horizontal branch stars). A third of the 
main sequence stars are consistent with having solar metallicity, thus 13% of 15**^ mag A 
stars are possible runaways. In our 2 Mq simulation, however, the fraction of runaways is 
six times larger than than the observed f raction in sta r coun ts at (7 = 15 (Figure [TOl) . Thus, 
if the main sequence A stars observed by iBrown et al. J2008h are all runaways, the absolute 
fraction of A stars ejected as runaways i s no larger than 0.13/6 ~ 2%, in agreement with the 
upper limit estimated by IStetsonI (Il98ll ). 



5.3. Runaways in the Inner Halo 



We now broaden our discussion and look at the distribution of runaways in the context 
of Galactic structure. 



Spatially, 80% of runaways in our simulations are located at |6| < 15° (Figure [TOj) . a 
population an observer migh t call "thick disk." The Galactic thick disk has an observ ed scale 
height of 0.75 - 1.0 kpc (i.e. Isiegel et al.ll2002l : iGirard et allbood : ljuric et al.l 120081 ) . and it 
dominates the number density of stars in the region 1 < |z| < 5 kpc. We note that claims 
of unusually large th ick disk scale heights are possibly confused with the inner stellar halo 
(IKinman et al.ll2009l ). By comparison, the number density of runaways in our simulations, 
selected with 7 < rcyi < 9 kpc and 1 < |2;| < 5 kpc, is well fit by an exponential distribution 
with a vertical scale height = 0.54 ± 0.01 kpc. 

Kinematically, runaways are a hot, rotating population analogous to the thick disk and 
inner stellar halo. In the region away from the plane |6| > 15°, simulated runaways have a 
mean 135 km s~^ component of v elocity in the directio n of Galactic rotation, comparable to 
that observed for the thick disk (jChiba fc Beersll2000l ). Runaways also have a 130 km s~^ 
heliocentric radial velocity dispersion, essentially identical to the velocity dispersion of inner 
halo stars with the same apparent magn i tude. The Milky Way inner halo, as described by 
CaroUo et al.l (120071 ) and iMorrison et al.l (120091 ) . has a small ~25 km s~^ prograde rotation 
and a ID velocity dispersion of ~120 km s~^. Thus runaways have similar kinematics to the 
Galactic thick disk and inner halo and are difficult to identify by radial velocity alone. 
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The m e an m etallicity of the inner halo is [Fe/H]= —1.6 (jCaroUo et al.l 120071 ). but 
Ivezic et al.l (120081 ) report solar metallicity stars up to 5 kpc above the plane. Perhaps run- 
aways can contribute to the high metallicity population and to the small prograde rotation 
observed in the inner stellar halo. 

The inner halo is much too metal poor to be composed entirely of runaways, yet metal 
rich, short-lived runaways should be present in the halo of the Milky Way. The deaths of 
massive runaways must necessarily result in metal enrichment and energy input into the halo. 
This conclusion applies for all star-forming galaxies, at all redshifts. In particular, runaways 
should be more abundant early in the evolution of a galaxy when the star formation rate 
is larger. Thus, the distribution of runaways may have important implications for feedback 
processes. 



6. Conclusions 

Runaway stars are an interesting class of objects because they connect star formation 
in the disk with the h alo of the Milky Way. We explore these connections by using the 



Portegies ZwartI (120001 ) distribution of binary-supernova ejections to inject stars into the 
Galactic potential. We track the progress of these stars from the Galactic disk to the Galactic 
halo and derive simulated catalogs of runaways. 

We show that the velocity and spatial distributions of runaways depend on the Galactic 
potential and stellar lifetime. All runaways have a flattened spatial distribution, with higher 
velocity runaways at lower Galactic latitudes. Massive runaways do not live long enough to 
reach the outer halo. Thus, massive runaways are more concentrated towards the Galactic 
center and in the plane of the disk than low mass runaways. 

Kinematically, runaways are a hot, rotating population of stars with dynamical proper- 
ties between the thick disk and the halo. In the solar neighborhood, runaways with masses 
of 1.5-3 Mq have scale heights and rotation velocities comparable to the thick disk, and ve- 
locity dispersions comparable to the inner stellar halo. Although they do not have a unique 
signature in radial velocity surveys, runaways are overwhelmingly located at low Galactic 
latitudes and at bright apparent magnitudes. Our results suggest an upper limit of 2% of 
A-type stars ejected as runaways. 

The kinematics of the unbound runaway HD 271791 provides an interesting comparison 
with the simulations. For runaways that reach distances of 60-100 kpc from the Galactic 
Center, the simulation s predict a peak in the initial distance of bound (unbound) stars at 



7-12 kpc (10-15 kpc). iHeber et al.l (120081 ) conclude that HD 271791 originated from 12 kpc 
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< r < 16 kpc (iHeber et al.ll2008l ). Clear tests of the simulations require larger samples of 
runaways with high quality proper motions and good estimates for their starting locations 
in the disk. 

Radial velocity surveys for high velocity outliers are unlikely to confuse runaway stars 
with HVSs. Among ejected stars with velocities exceeding 400 km s~^, runaways (1) are 
brighter {g < 18), (2) have smaller velocities {v < 450 km s~^), and (3) are more concentrated 
to lower galactic latitudes (b < 35°) than HVS with typical g > 18, v ^ 400-700 km s~^, and 
random b in the range 30° < b < 90°. We estimate that at most ~ 1 runaway contaminates 
the Brown et al. sample of HVSs. 

Future progress on the theory of runaways requires predictions for the velocity distribu- 
tion from the dynamical ejection mechanism. Although this process can yield higher ejection 
velocities than the binary-supernova ejection mechanism, uncertainties in the stellar merger 
rate during close encounters complicates calculations of a realistic maximum ejection velocity 
and a distribution of ejection velocities. 

For both mechanisms, better estimates of predicted rates are nee ded to constrain predic- 



tions for the frequency and kinematics of runaways. Extending the iPortegies ZwartI (120001 ) 
simulations to lower mass stars would improve our estimate for the rate of runaways from 
the binary-supernova ejection mechanism among A-type stars. Simulations of ensembles of 
dense star clusters would yield ejections rates for the dynamical ejection mechanism. 

Measurements of proper motions and metallicities of inner halo stars can also improve 
our understanding of runaways. Identifying the fraction of halo stars with roughly solar 
metallicity can yield a better estimate of the frequency of runaways. Proper motions and 
radial velocities place clear constraints on the origin of runaways for comparison with theo- 
retical simulations. 
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Table 1. Median Radial Velocities of Runaway Stars 


a 


bf (deg) 


Galactocentric distance r 
< 20 kpc 20-40 kpc 40-60 kpc 60-80 kpc 80-100 kpc 


100-120 kpc 




= 1.5 Mq, tms = 2.9 Gyr, rinit = 10 kpc, Vej= 400 km s"^ 


0-30 
30-50 
50-90 


20 15 26 25 31 
69 -35 -46 -37 10 
-13 94 92 51 34 


44 
23 
23 


= 2 M0, tms = 1-2 Gyr, ri„it = 10 kpc, Vej= 400 km s'^ 


0-30 
30-50 
50-90 


45 54 83 78 150 
110 44 -33 -27 40 
10 93 111 62 35 


170 
58 
24 


= 3 M0, tms = 350 Myr, ri„it = 3 kpc, Vej= 400 km s'^ 


0-30 

30-50 
50-90 


86 122 142 133 150 

-36 22 40 37 
100 33 


158 


M* = 3 Mq, tms = 350 Myr, ri„it = 10 kpc, Vej= 400 km 


0-30 

30 50 
50-90 


112 190 240 275 295 

222 266 230 210 205 
166 190 152 115 130 


304 
230 


= 3 Mq, tms = 350 Myr, li^it = 30 kpc, Vej= 400 km 


0-30 

30-50 
50-90 


-36 112 314 339 363 
-26 165 281 321 330 
60 145 226 254 259 


384 
333 
258 


= 3 Mq, tms = 350 Myr, ri„it = 10 kpc, Vej= 350 km 


0-30 

30-50 
50-90 


90 155 207 228 232 
197 191 104 125 155 
98 127 74 42 


243 
176 
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Table 1 — Continued 



Galactocentric distance r 
bf (deg) < 20 kpc 20-40 kpc 40-60 kpc 60-80 kpc 80-100 kpc 100-120 kpc 







4 Mq, tms 


= 160 Myr, rinit = 


10 kpc, 


, Vcj= 400 km 


30 


223 


320 


361 


382 


404 


30 50 


249 


308 


295 


319 




50-90 


201 


232 


221 








= 


5 Mq, tms 


= 95 Myr, ri„it = 


10 kpc, 


Vej= 400 km s-i 


0-30 


254 


388 


418 






30-50 


260 


335 


359 






50-90 


217 


270 


301 








= 


6 Mq, tms 


= 65 Myr, rinit = 


10 kpc. 


Vej= 400 km 


0-30 


271 


418 








30-50 


268 


351 








50-90 


239 


292 









^Results for fiducial simulations of 10^ stars with one initial starting radius, stellar mass, 
and ejection velocity for all stars. Each entry lists the median Vrad, measured in the 
Galactocentric reference frame, as a function of r and bf for stars with the listed initial 
conditions. For entries without data, simulations with 10^ trials do not produce stars with 
these combinations of r and bf. 
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Table 2. Radial Velocity Dispersions of Runaway Stars ^ 









Galactocentric distance r 






bf (deg) 


< 20 kpc 


20-40 kpc 


40-60 kpc 60-80 kpc 80-100 kpc 


100-120 kpc 








iVl^ — ±.0 iVl0 


5 Iinit — J^pt^? Vej — ^UU Kill o 


-1 




0-30 


193 


204 


177 146 


140 


144 


30-50 


200 


212 


190 168 


154 


135 


50-90 


184 


158 


121 85 


70 


46 








T*. — in If-r*/-" T7- . — Ann L-m a~ 
■^init — ■'^P^; ^ej — rtKJKJ Kill o 


-1 




0-30 


192 


199 


172 172 


174 


154 


30-50 


190 


220 


200 162 


134 


110 


50-90 


187 


172 


133 90 


69 


46 






AT _ O AT„ 
iVi* — O iVl0, 


V. .. — L-T^r* ■\T ■ — /inn If-m a 
■•■init — " -^P^? *ej — ^yjyj Kill o 


1 




0-30 


179 


171 


108 64 


30 


4 


30-50 


206 


134 


72 27 






50-90 


140 


78 












= 3 Mq, 


Tinit = 10 kpc, Vej= 400 km s" 


-1 




0-30 


180 


201 


157 117 


83 


54 


30 50 


155 


165 


123 84 


55 


33 


50-90 


160 


136 


87 48 


23 








= 3 Mq, 


Tinit = 30 kpc, Vej= 400 km s" 


-1 




0-30 


257 


180 


137 128 


109 


89 


30-50 


179 


120 


104 100 


82 


61 


50-90 


101 


87 


94 77 


54 


34 






= 3 Mq, 


Tinit = 10 kpc, Vej= 350 km s" 


-1 




0-30 


165 


183 


134 92 


58 


34 


30-50 


143 


148 


101 62 


32 


8 


50-90 


138 


112 


61 21 
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Tabic 2 — Contimicd 









Galactocentric distance r 


bf (deg) 


< 20 kpc 


20-40 kpc 


40-60 kpc 60-80 kpc 80-100 kpc 100-120 kpc 








= 4 Mq, 


Tinit = 10 kpc, Vej= 400 km s ^ 


0-30 


148 


141 


85 41 8 


30-50 


124 


108 


58 27 


50-90 


111 


83 


36 






= 5 Mq, 


Tinit = 10 kpc, Vej= 400 km 


0-30 


130 


103 


46 


30-50 


100 


74 


31 


50-90 


82 


52 


10 






= 6 Mq, 


Tinit = 10 kpc, Vej= 400 km 


0-30 


127 


85 




30-50 


88 


61 




50-90 


68 


35 





Results for fiducial simulations with one initial starting radius and ejection velocity for 
each stellar mass. Each entry lists the ID radial velocity dispersion as a function of r and 
bf for stars with the listed initial conditions. For entries without data, the simulations do 
not produce stars with these combinations of r and bf. 
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Table 3. Median Radial Velocities of Runaway Stars Ejected from the Galactic Disk' 



Galactocentric distance r 
Stellar Mass < 20 kpc 20-40 kpc 40-60 kpc 60-80 kpc 80-100 kpc 100-120 kpc 



h < 30° 


1.5 Mq 


1 


3 


7 


9 


15 


16 


2 Mq 


2 


8 


17 


29 


28 


42 


3 Mq 


8 


23 


84 


136 


188 


239 


4 Mq 


20 


114 


239 


332 


411 




5 Mq 


53 


202 


361 








6 Mq 


72 


246 


355 








30° < 6 < 50° 


1.5 Mq 


4 


10 


14 


15 


16 


23 


2 Mq 


13 


17 


16 


22 


27 


38 


3 Mq 


50 


46 


73 


121 


172 


224 


4 Mq 


70 


115 


217 


309 






5 Mq 


103 


199 


334 








6 Mq 


144 


260 










b > 50° 


1.5 Mq 


22 


48 


54 


52 


68 


49 


2 Mq 


37 


64 


70 


71 


65 


48 


3 Mq 


72 


88 


88 


108 


152 


208 


4 Mq 


97 


127 


194 


285 






5 Mq 


128 


187 


287 








6 Mq 


169 


245 











^Results for 10^ stars ejected with velocity Vinit from r = rinit in the Galactic Disk. The 
ejection velocity and position are chosen from probability distributions (Eq. (SHHl) described in 
the main text. Each column lists the median Vj-ad for 1.5-6 Mq stars with r in the listed range. 
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Table 4. Radial Velocity Dispersions of Runaway Stars Ejected from the Galactic Disk' 









Galactocentric distance 


r 




Stellar Mass 


< 20 kpc 


20-40 kpc 


40-60 kpc 


60-80 kpc 


80-100 kpc 


100-120 kpc 




b < 30° 


1.5 Mq 


102 


127 


129 


124 


116 


115 


2 Mq 


99 


125 


126 


124 


122 


112 


3 Mq 


94 


125 


109 


86 


66 


50 


4 Mq 


93 


102 


71 


44 


32 




5 Mq 


94 


102 


53 








6 Mq 


104 


125 


62 








30° < & < 50° 


1.5 Mq 


121 


127 


117 


116 


117 


111 


2 Mq 


120 


127 


125 


118 


106 


97 


3 Mq 


118 


118 


96 


73 


55 


40 


4 Mq 


102 


90 


58 


34 






5 Mq 


88 


74 


41 








6 Mq 


83 


67 










b > 50° 


1.5 Mq 


121 


114 


96 


84 


87 


82 


2 Mq 


121 


110 


99 


91 


88 


80 


3 Mq 


117 


104 


79 


56 


40 


25 


4 Mq 


103 


76 


42 


15 






5 Mq 


82 


54 


16 








6 Mq 


67 


36 











As in Table [3] for the radial velocity dispersion. 
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Table 5. Vertical Velocity Dispersions of Runaway Stars Ejected from the Galactic Disk' 



Stellar Mass 


< 20 kpc 


20-40 kpc 


Galactocentric distance 
40-60 kpc 60-80 kpc 


r 

80-100 kpc 


100-120 kpc 


1.5 Mq 


65 


60 


49 


42 


36 


35 


2 Mq 


64 


61 


51 


47 


42 


37 


3 Mq 


67 


67 


50 


41 


41 


42 


4 Mq 


70 


62 


57 


58 


52 




5 Mq 


69 


69 


68 








6 Mq 


71 


78 


74 









As in Table [3] for the velocity dispersion perpendicular to the disk. 
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Table 6. Vertical Scale Height (in kpc) of Runaway Stars Ejected from the Galactic Disk' 



Stellar Mass 


< 20 kpc 


20-40 kpc 


Galactocentric distance 
40-60 kpc 60-80 kpc 


r 

80-100 kpc 


100-120 kpc 


1.5 Mq 


1.5 


4 


9 


15 


20 


25 


2 Mq 


1.5 


4 


10 


15 


25 


30 


3 Mq 


1.5 


4 


10 


20 


30 


40 


4 Mq 


1.5 


4 


15 


25 






5 Mq 


2 


4 


15 








6 Mq 


2 


4 


15 









As in Table [3] for the vertical scale height. 
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Table 7. Median Properties for Unbound Runaway Stars. 



Stellar Mass 


r (kpc) 


Vrad (km S ^) 


bf (deg) 


1.5 Mq 


355 


210 


9 


2 M0 


210 


263 


9 


3 M0 


82 


360 


10 


4 Mq 


43 


400 


9 


5 Mq 


30 


410 


8 


6 Mq 


23 


410 


7 
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r (kpc) 



Fig. 1. — Scatter plots of simulated runaways ejected with rinit = 10 kpc, Vej = 400 km s~ , 
and m = 3 Mq. Each point is a (r,Vi.ad) pah derived from a Monte Carlo simulation of 10^ 
stars. Upper panel: points ejected into the Galactic disk (6j < 30°) are color-coded according 
to the ejection angle 6i relative to the direction of Galactic rotation. Blue (yellow) points 
indicate stars ejected with 9i = —tx/A — > 7r/4 {6i = Sn/A 57r/4). Cyan and magenta 
points indicate stars ejected at intermediate angles, 9i = ±(7r/4 7r/2) for cyan points 
and 9i = ±(7r/2 37r/4) for magenta points. Stars ejected along (against) the direction of 
Galactic rotation reach the largest (smallest) distances in the Galactic halo. Lower panel: 
points are color-coded according to bi, the ejection angle relative to the Galactic plane, with 
bi < 30° in cyan, 30° < 6, < 50° in magenta, and bi > 50° in yellow. Stars ejected into the 
plane (halo) reach the largest (smallest) distances from the Galactic center. 
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r (kpc) 

Fig. 2. — Scatter plots of simulated runaways ejected with fixed initial radius rinit, ejection 
velocity Vgj, and stellar mass m. Each point is a (r,Vrad) pair derived from a Monte Carlo 
simulation of 10^ stars. In each panel, the colors of the points indicate the Galactic latitude 
at the end of the simulation; cyan: b < 30°, magenta: 30° < b < 50°, and yellow: b > 50°. 
Upper left panel: results of the fiducial simulation with rjnit = 10 kpc, Vej= 400 km s~^, and 
m = 3 Mq. Runaways ejected parallel to the plane of the Galactic disk achieve higher Vrad 
at larger r than runaways ejected perpendicular to the disk plane. Upper right panel: results 
for rinit = 30 kpc. Runaways ejected at larger rinit achieve larger Vrad at larger r. Lower left 
panel: results for Vej= 350 km s~^. Runaways ejected with smaller Vgj have smaller Vrad at 
all r. Lower right panel: results for m = 4 Mq. Smaller main sequence lifetimes do not allow 
more massive stars to reach large r. 
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50 
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Fig. 3. — Comparison of (r, Vrad) diagrams for halo stars, HVSs, and runaways. Stars are 
color-coded according to their Galactic latitude; cyan: &/ < 30°, magenta: 30° < &/ < 50°, 
and yellow: 6/ > 50°. The top panels shows results for 3 Mq (left panel) and 4 Mq (right 
panel) runaways ejected from the Galactic disk with probability distributions for Vej (Eq. E]) 
and for rinit (Eq. [5]). Longer- lived lower mass runaway stars reach larger r than shorter-lived 
more massive runaways. Runaways ejected into the disk achieve larger r than stars ejected 
perpendicular to the disk. The lower left panel plots the expected velocity distribution of 
halo stars with a velocity dispersion of 110 km s^^ and a oc density law. Although 
nearby halo stars can have large Vrad, distant runaways have much larger Vrad than distant 
halo stars. The lower right panel shows the expected velocity distribution of HVSs ejected 
from the Galactic Center. At each r, the lower limit to Vj-ad for HVSs follows the predicted 
lower limit for runaway stars. However, at all r and bf, HVS have much larger v^ad than 
runaway stars. 
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Fig. 4. — Predicted surface density distributions as a function of cylindrical radius {riyi = 
x'j + y'j) for runaways ejected with a range of velocities (Eq. [6]) from an exponential disk (Eq. 
[5]). The dot-dashed black line plots the initial density distribution. The colored lines plot 
surface density distributions for stars of different masses, as indicated in the legend. For all 
stars, the density rises slowly from 1-3 kpc, follows a power law (S oc r~^'^) at intermediate 
radii, and then falls exponentially at large radii. 
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z (kpc) 

Fig. 5. — Predicted natural log of the relative number density as a function of 2;, the height 
above the disk plane, for 2 Mq runaway stars. The curves plot the number density in six 
cylindrical radius bins, as indicated. In each bin, the density distribution has a central core 
with a scale height of 300-1000 pc and an extended halo with a scale height of 2-30 kpc. 
The z scale height of the extended halo increases with increasing cylindrical radius. 
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Fig. 6. — Cumulative probability distributions as a function of initial disk radius for 2 Mq 
runaway stars. Among runaways that reach final radii r > 60 kpc (dashed line), r > 80 kpc 
(solid line), and r > 100 kpc (dot-dashed line), roughly 50% are ejected from inside the solar 
circle. Among unbound runaways (double dot-dashed line), roughly 50% are ejected from 
rinit ^ 12 kpc. 



10 20 30 40 

b„ (degrees) 



Fig. 7. — Cumulative probability distributions as a function of final galactic latitude for 
unbound runaway stars. The legend indicates the stellar mass for each curve. The median 
galactic latitude ranges from bf^rned = 7° for 6 Mq stars to bf med = 9°-10° for 1.5-3 Mq 
stars. All unbound runaways have bf < 30 °. 
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r (kpc) 

Fig. 8. — Cumulative probability distributions as a function of distance for unbound runaway 
stars. The legend indicates the stellar mass for each curve. Lower mass runaways reach larger 
r than more massive runaway stars. 
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Fig. 9. — Predicted distributions of runaways and HVSs in a heliocentric reference frame. 
The top panels plot number fractions of simulated runaways {solid lines) and HVSs {dashed 
lines) as a function of Galactic longitude, latitude, and apparent magnitude. Bottom panels 
plot the 50**^, 90**^, and 99*'' percentile heliocentric radial velocities of simulated runaways 
and HVSs. The spatial and velocity distributions of 3 Mq runaways and HVSs are very 
different: runaways are apparently bright and concentrated to the Galactic plane; HVSs are 
apparently faint and common in the halo. 
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Fig. 10. — Comparison of number fractions of simulated 1.5, 2 and 3 runaways in 
the SDSS survey region {blue lines) with SDSS star counts {red lines). Relative to the 
Galactic stellar population, runaways have the highest contrast at low latitudes and at 
bright magnitudes. 



